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Abstract 

I calculate the first correction to the thermal distribution function of an expanding gas due to 

shear viscosity. With this modified distribution function I estimate viscous corrections to spectra, 
Ch ; 

elliptic flow, and HBT radii in hydrodynamic simulations of heavy ion collisions using the blast 
wave model. For reasonable values of the shear viscosity, viscous corrections become of order one 
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when the transverse momentum of the particle is larger than 1.7 GeV. This places a bound on 
the pt range accessible to hydrodynamics for this observable. Shear corrections to elliptic flow 
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cause V2{pt) to veer below the ideal results for px ~ 0.9 GeV. Shear corrections to the longitudinal 
HBT radius R\ are large and negative. The reduction of R\ can be traced to the reduction of 
the longitudinal pressure. Viscous corrections cause the longitudinal radius to deviate from the 
scaling which is observed in the data and which is predicted by ideal hydrodynamics. The 
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correction to the sideward radius Rg is small. The correction to the outward radius Rq is also 
negative and tends to make Ro/Rs ~ 1- 
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I. INTRODUCTION 



One of the most exciting results of the Relativistic Heavy Ion Collider (RHIC) is the 
observation of collective motion. In particular, the experiments have measured a large 
elliptic flow in non-central collisions [3, Q, 0, U, 0] . Elliptic flow is quantified with the second 
harmonic of the azimuthal distribution of produced particles 



j\d<P cos(20) - v 

-7v dy p t dpt d(f> 



v 2 ( PT ) = (cos(20)) PT = - - ) (1) 



where <fi is the measured relative to the reaction plane, ^(pt) rises strongly as a function 
of transverse momentum up to ~ 1.5 GeV. One interpretation of the observed flow is 
that hydrodynamic pressure is built up from the rescattering of produced secondaries and 
pressure gradients subsequently drive collective motion. A strong hydrodynamic response 
is possible if the sound attenuation length T s = f^;, is significantly smaller than the 
expansion rate, ~ r. (In the formula T s = fj^j, V is the shear viscosity, e the energy 
density and p the pressure.) Estimates based upon perturbation theory give T s ~ r and 
indeed thirty times the perturbative 2-2 cross sections are needed to obtain the observed 

n 

elliptic flow [f|. However, these perturbative estimates are uncertain. In an example of a 
strongly coupled gauge theory where calculations are possible (N=4 SUSY YM), T s is in 
fact approximately 2-4 times smaller compared to perturbation theory ?| (see also Section 

Ideal hydrodynamics (r s = 0) has been used to simulate heavy ion reactions and readily 
reproduced the observed elliptic flow and its dependence on centrality, mass, beam energy 
and transverse momentum 0, E| . However ideal hydrodynamics failed in several respects. 
First, above px ~ 1.5 GeV the observed elliptic flow does not increase further as predicted 
by hydrodynamics. Additionally, the single particle spectra deviate from hydrodynamic 
predictions above p? k, 1.5 GeV. Second, the observed HBT radii are significantly smaller 
than predicted by ideal hydrodynamics ^0, 11, In particular, the longitudinal radius Rl 
is 50% smaller than the ideal hydrodynamic result. Further, the ratio between the outward 
(Ro) an d sideward (Rs) radii is observed to be approximately one while ideal hydrodynamics 
predicts Ro/Rs » 1.3 (ifj. 

The domain of applicability of hydrodynamics can be answered quantitatively by calcu- 
lating the first viscous correction to ideal hydrodynamic results. The effect of viscosity is 
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twofold. First, viscosity changes the solution to the equations of motion. Second, viscosity 
changes the local thermal distribution function. This effect was first investigated in heavy 
ion physics by Dumitru 13J. The purpose of this work is to consider the effect of a modi- 
fied thermal distribution function on spectra, elliptic flow, and HBT radii. Thus this work 
delineates the boundaries of the hydrodynamic description as applied to relativistic heavy 
ion collisions. 



II. VISCOUS CORRECTIONS TO A BOOST INVARIANT EXPANSION 

First consider a baryon free viscous boost invariant expansion with a vanishing bulk 
viscosity, but a non-zero shear viscosity, 77. Note throughout this work we denote the space- 
time rapidity as r] s and the viscosity as rj. Unlike for ideal hydrodynamics where entropy is 



conserved, the entropy per unit space-time rapidity ts increases as a function of r = \Jt 2 — z 2 

HQ, HQ 

d ( TS ) _ iv (9 s 

dr ~tT- {Z) 

For hydrodynamics to be valid, the entropy produced over the time scale of the expansion 
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t (to wit, t- 2 ^) must be small compared to the the total entropy, (ts). This leads to the 
requirement that 

- « 1 , (3) 
r 

where we have defined the sound attenuation length 

r s = % ■ (4) 

sT v ; 

T s is approximately the mean free path and therefore the condition T s /t <d 1 is just the 
statement that the mean free path be small compared to the system size. The name "sound 
attenuation length" follows from the dispersion relation for a sound pulse u = c s 
where c 2 s = (||) is the squared speed of sound. In the remainder of this section, I gather 
estimates for T s in the Quark Gluon Plasma (QGP). For similar estimates in the hadron gas 
see Q. 

The shear viscosity has been determined in the perturbative QGP only to leading log 



accuracy 



ial2o|. To leading log (g : ) the shear viscosity with two light flavor is given by 
86.473 pr iogfa*-!) • With the entropy of the QGP, s = 37^-T 3 and setting a s — ► \ and 
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log(g x ) — > 1 the sound attenuation length in perturbation theory is 



(5) 



r s \ 1 

-1 = 0.18— 

T J Pert. TT 

Estimates of evolution time scales give tT ~ 1. The value of T s /t is sensitive to the value 
of a s . 

This perturbative estimate of T s is clearly uncertain and assumes that a s m 1/2 and 
that log(g _1 ) is a large number. Recently the shear viscosity was evaluated in a strongly 
coupled gauge theory, N = A SUSY YM using the AdS/CFT correspondence 1711 . The shear 



viscosity is given by rj = ^N 2 T 3 [?| and the entropy is given by s = ^N^T 3 |2l|. Thus in 
this strongly coupled field theory T s is 

T / AdS/CFT 3-KTT 

which is 2-4 times smaller than the corresponding perturbative estimate depending. 

Finally, I compare these theoretical estimates of T s to the value abstracted from Monte 
Carlo simulations of RHIC collisions performed by Gyulassy and Molnar (GM) Q|. GM 
modeled the heavy ion reaction as a gas of massless classical particles suffering only 2 — > 2 
elastic collisions with a constant cross section in the c.m.s frame, ^ = f 2 -. When particle 
number is conserved, T s is given by a more complicated formula which reflects the coupling 
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between the energy and number densities 
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e + p 



or 



e + p e + p \dT J n 
where k is the thermal conductivity. For the GM gas, c 2 s = 
to as before. The shear viscosity in the GM gas is r] 
directly proportional to the mean free path 



\ dT Jn c s \9n 



(7) 
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0.421- 
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nT and Y s reduces 
23( . Therefore T s is 

(8) 



In order to achieve a reasonable agreement with the measured elliptic flow, GM required a 
transport opacity of x ~ 20-^40. This transport opacity was reached when the cross section 
was o"o ~ 10 + 20 mb and the number of particles was ^ ~ 1000 at proper time r = 0.1 fm. 
The initial density of particles is n = /(t ttR 2 ). Substituting R ps 5.5 fm we obtain 

T ' GM 



0.02^-0.04 . 



(9) 



This is smaller by a factor of three or more than even the AdS/CFT estimate assuming that 
tT ~ 1. The physical mechanism for such a small viscosity remains unclear. 

The sound attenuation length is uncertain. In what follows we take — = \ and calculate 

° TO 

viscous corrections to the observed spectra, elliptic flow, and HBT radii. In summary, 
perturbation theory finds T s /r pa 0.18, strongly coupled supersymmetric field theory finds 
T s /r pa 0.11, and phenomenology finds T s /t pa 0.03. 



III. VISCOUS CORRECTIONS TO THE DISTRIBUTION FUNCTION 



Viscosity modifies the thermal distribution function. The formal procedure for deter- 
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re viscous corrections to the thermal distribution function is given in the refer- 



24j . In general, for a multi-component gas the viscous correction is different for 



each component. For simplicity, we will consider a single component gas of "pions" with 
= 140 MeV. The basic form of the viscous correction can be intuited without calculation. 
First write f(p) = f Q + 5f, where f {^) = p-J/t^ is the equilibrium thermal distribution 
function and Sf is the first viscous correction. Sf is linearly proportional to the spatial gradi- 
ents in the system. Spatial gradients which have no time derivatives in the rest frame and are 
therefore formed with the differential operator V M = {g^ v — u p u v )d v . For a baryon free fluid, 
these gradients are V Q T, V Q w a , and (VaUp), where (V a up) = V a up + V pu a — |A Q , j aV 7 M 7 . 
V a T can be converted into spatial derivatives V ' a up using the ideal equations of motion and 
the condition that T^ u u v = eu M j^J. V a u a leads ultimately to a bulk viscosity and will be 
neglected in what follows. Finally, (V a up) leads to a shear viscosity. If Sf/f Q is restricted 
to be a polynomial of degree less than three in p^, then the functional form of the viscous 
correction is completely determined, 

/ = / (i + ^V(v Q ^)). (io) 

For a Boltzmann gas this is the form of the viscous correction adopted in this work. The 
factor of 2 in ^ is inserted for later convenience. For Bose and Fermi gasses the ideal 
distribution function in Eq. El is replaced with f {l ±f ) The correction described 

here is precisely the "first approximation" of reference 24| and the "one parameter ansatz" 
for a variational solution of reference 19]. The "one parameter ansatz" reproduces the full 
result to the 15% level. 
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The coefficient C in Eq. El can be reexpressed in terms of the sound attenuation length. 
Indeed, substituting / to determine the stress energy tensor 



we find 

d 3 p 



J (2^ W//o(1 + /o) 



(V a u p ) . (12) 



2T 3 

The quantity in square brackets is a fourth rank symmetric tensor and consequently can be 
written in terms of A^ v = g^ v — u^u u and u^. Thus, 

2T 3 " / 7^^^PV/ (l + / )^^ (13) 

+a 2 (A^A aP + A^ a A u/3 + A^A™) . 

Substituting Eq.[T3] into Eq.[T2] and using the identities u a (V a Up) = (V a up) = 
A al3 {V a up) = 0, we find 2ct 2 = r\. To determine the coefficient a 2 , contract both sides 
of Eq.IUwith 

and evaluate the resulting expression in the local rest frame. The result for the viscosity is 

1 = Mr f T^L/oO + /•) IpI 4 • (15) 



90 T 3 J (27r) 3 £- 

For a Boltzmann gas / G (l + / G ) is be replaced with / (^r) = e~^~ and the integrals can 
be performed analytically. Comparing the resulting expression to the entropy of an ideal 
Boltzmann gas (see e.g. [2 5] ) we find C = 2 . For a massless Bose gas the integrals can again 
be performed analytically and C = 90 ^ 5 ) 7 ~ 1-04^. For a massive Bose gas, the integral 
was performed numerically and C varies monotonously between these two limiting cases. 
Therefore up to a few percent, we have C = -, and the viscous correction 5f is 

*/ = ^/»(i + /o)pVM • 

IV. VISCOUS CORRECTIONS TO A BJORKEN EXPANSION 

Before considering the viscous corrections to more general hydrodynamic expansions, let 
us consider a simple Bjorken expansion of infinitely large nuclei without transverse flow. At 
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mid space-time rapidity the stress energy tensor is at time r a is given by [1 
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(16) 



p— g^-, while the transverse pressure is increased 



p 


\o 

where, denotes the ideal stress energy tensor diag(e, p,p,p), Thus, the longitudinal 
pressure is reduced by the expansion, T 
by the expansion, T x 

The difference between the longitudinal and transverse pressures is reflected in the pr 
spectrum of thermal distribution. Since the transverse pressure (T xx ) is increased by |^-, the 
particles are pushed out to larger p^. Armed with the modified thermal distribution function, 
the Cooper Frye formula 2(| gives the thermal spectrum of particles in the transverse plane 
at proper time r D 

d 2 N 1 



d 2 N^ 



+ 



d?pT dy 
d 2 N^ 



(2tt) 3 
1 



p^fo + Sf 



(17a) 
(17b) 



d 2 pTdy d 2 pxdy (2tt) 3 

Here dT,^ is the oriented space-time volume. Substituting into Eq. [T7|(see Appendix |B} we 
obtain the the ratio between the viscous correction (5 dN = -jjjr—}f ) and the ideal spectrum 

{dN (0) - diV(0) N 



cPpTdy i 



5dN 



4r \\t) \T J 2 \K!(^ 



k 3 (t: 



- 1 



Using the asymptotic expansion for the modified Bessel functions, we have for large trans- 
verse momenta, 



5dN 



4:T n 



Pry 
T ) 



(18) 



As promised, the larger transverse pressure drives pushes the corrected spectrum out to 
higher transverse momenta. For a Bjorken expansion without transverse flow, this formula 
also indicates at what transverse momentum the hydrodynamic description of px spectra is 
applicable. For y- ~ |, and T = 200 MeV the ratio between the ideal spectrum and the 
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FIG. 1: (a) Thep z distribution of particles with coordinate-space rapidity rj s = 0, with and without 
viscous corrections, (b) The z distribution of particles with momentum-space rapidity y = 0, with 
and without viscous corrections. The curves are drawn for a Bjorken expansion without transverse 
flow at t = 7fm for a Boltzmann gas with temperature, T = 160 MeV, m = 140 MeV. The 
transverse momentum is fixed, pt = 400 MeV. The viscous correction is linearly proportional to 

r s /r Q . 

correction becomes of order one for p™ ax m 800 MeV. We shall see in the next section that 
this upper bound on the domain of hydrodynamics is significantly larger p™ ax « 1.5 GeV 
once the transverse expansion is included in the flow profile. 

We have already noted that the longitudinal pressure is reduced by the expansion, T zz = 
p — The reduction in the longitudinal pressure is ultimately responsible for a reduction 
in the longitudinal radius measured by Hanbury-Brown Twiss interferometry. Since the 
longitudinal pressure is reduced due to the expansion, the distribution in p z at mid space- 
time rapidity (r] s = 0) is narrower. This is illustrated in Fig. HJa) for a fixed transverse 
momentum px = 400 MeV. Due to boost invariance the p z distribution at r] s = is directly 
related to the z distribution at y = Specifically, for fixed transverse momentum, " A 
is a function of \y — rj s \, which leads to the relation 



dydri s 



rriT 



dN 
dp z dr] 



ri=0 



dN 
' dydz 



(19) 
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It follows that the z distribution at mid momentum-space rapidity is narrower as indicated 
in Fig. HJb). The width of this z-distribution, is related to the longitudinal radius that is 
measured by HBT interferometry (see e.g. \^\). 

To understand this result analytically we must calculate the width of z distribution for 
a simple Bjorken expansion of a Boltzmann gas at proper time r Q . Let us quickly recall 
the definitions of the HBT radii. The source function S(x, K) for on shell pion emission is 
defined such that 

E k~^ = \ d 4 xS(x,K) (20) 



d 3 K 

where Ek = K° = a/K 2 + m 2 . Averages with respect to the source function are defined as 
(a) K = J d 4 xaS(x,K)/ J d 4 xS(x,K). To a good approximation (see e.g. Ref j^). certain 
spatial and temporal variances of the source function can be determined from the Bose- 
Einstein correlations between pion pairs at small relative momenta. For a boost invariant 
and rotationally invariant source, we can assume without loss of generality that the pair 
momentum points in the x direction (i.e. K = (K x ,K y ,K z ) = (Kt,0, 0)). Then the 
following variances can be determined from HBT measurements 

R 2 o(K T ) = {{x-v K tf) KT (21) 
Rl(K T ) = (f) KT (22) 
Rl(K T ) ee (~z 2 ) Kt , (23) 

where vk = Kt/Ek and for example x = x — (x). Comparing Eq. ^Jand Eq. w e see 
that in this work the source function is confined to a freezeout surface and therefore the 
averages are understood to mean 

_ S s K>dE f anx,K) 

The assumption of a sharp freezeout surface is clearly unrealistic. In general there is a 
transition region from hydrodynamics to the Knudsen limit. Within ideal hydrodynamics 
this transition region can not be determined. Within viscous hydrodynamics, viscous terms 
become large (~ 1/2) and signal the transition. 

Armed with these formula, the computation of R\ for a boost invariant expansion is 
straight forward. We have 



Rl(K t ) = {z ) Kt = 5K ^ f ^ K) ■ 



(25) 



Substituting / = f a + 5f, expanding to first order in 5f, and performing the integrals (see 



Appendix |Bj) we find the viscous correction 8R 2 L 



5Rl r s \6m T K 3 (^) _ , mr y 1 f K 3 (^) 



(i?l)(o) 



4 T K 2 (^) V T / 8 V^2(T) 



(26) 



where the (i?|)( ) is the ideal longitudinal radius [2£ 



(0) = 2 ^^M . (27) 

For the relevant range of T=f, the Bessel function expression in square brackets is large ^6—8. 
Accordingly, viscous corrections to the longitudinal radius are quite large (> 100%) and tend 
to reduce the radius relative to its ideal value. Including the transverse expansion reduces 
the viscous correction to 50% . Nevertheless, the viscous correction to the longitudinal 
radius remains large unless T s /r is significantly smaller than 0.1 . This formula and some 
caveats are discussed further in the next section. 

V. VISCOUS CORRECTIONS WITH TRANSVERSE EXPANSION 

To go further and illustrate the effect of viscosity on the observed spectra, elliptic flow 
and HBT radii of hydrodynamical models of the heavy ion collision, I generalize the blast 
wave model to include the viscous corrections of Eq.El The blast wave model provides of 
a simple parametrization of the flow of full ideal hydrodynamic simulations which assume 
boost invariance J^l 0] . The corrections described below are therefore indicative of similar 
corrections to these simulations. This is the reason for adopting the blast wave model here. 
The blast wave model also has been used to fit experimental data. The model provides 
a good description of spectra aad elliptic flow S B Q and provides a fair description of 
HBT radii for small Mr, Mr < 0.5 GeV [30]. However, for larger My the model does not 



reproduce the strong dependence on M T seen in the Ro and Rs radii 29, 31]. The blast 
wave model remains simply a model of the flow fields and ultimately a full viscous simulation 
is needed to estimate viscous effects. 

In the blast wave model of central collisions considered here, a hot pion gas is expanding 
in a boost invariant fashion and freezes out at a proper time r Q . In the transverse plane, 
the temperature is constant T Q = 160 MeV and the matter distribution is uniform up to a 



radius R Q . The transverse velocity rises linearly as a function of the radius, u r = u 

10 



o r 

Ro' 





Central (0-5%) 


Non-central (16-24%) 


T (MeV) 


160 


160 


Ro (fin) 


10 


7.5 


r Q (fm) 


7.0 


5.25 


u D 


0.55 


0.55 


u 2 





0.1 



TABLE I: Table of parameters used in the blast wave model described in the text. 
Summarizing, the hydrodynamic fields (T and m m ) are parameterized as 



T(T ,r] s ,r,, 
u r (r ,r] s ,r, , 



u 



T o e(R -r) 

T 

u — Q(R - r) 







(28a) 
(28b) 
(28c) 
(28d) 
(28e) 



The blast wave parameters are adjusted so that model with the ideal thermal distribution 
can approximately reproduce the spectra and HBT radii. Similar blast wave model fits 
have appeared ubiquitously in the heavy ion literature (see e.g. [29]). Then with the model 
parameters fixed, the viscous correction is calculated and compared to the ideal results. The 
model parameters for central collisions are recorded in Table|I] 

With the hydrodynamic fields specified, the viscous tensor (\7 a u 13 ') can be computed in 
a simple but lengthy calculation which is worked out in Appendix [XJ One technical point 
should be noted. In the viscous tensor {[V Q «^) time derivatives of the velocity appear. These 
time derivatives are converted into spatial derivatives using the ideal equations of motion 
which are sufficient to leading order in the viscosity. 

The spectrum of particles emerging from the freezeout oriented 3-volume is calculated 
by employing the Cooper-Frye formula, Eq. El These integrals are performed numerically 
in a straightforward fashion. Again relevant details are relegated to Appendix El The ideal 
spectrum of this blast wave model is typical of blast and is in rough agreement with pion 
data at RHIC. (See e.g. j^l for fits to data of this type.) In Fig. El the solid line shows 
the ratio of the viscous correction to the ideal spectrum. The dashed line shows the Bjorken 
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FIG. 2: The solid line shows the ratio between the viscous correction (5 dN = J^J^ ) and the 
ideal spectrum (dN^ = ) • The dashed line shows the Bjorken result without transverse flow 

given in Ea. ll81 The band indicates where the hydrodynamic description of the pt spectrum in the 
blast wave model can not be reliably calculated. The viscous correction is linearly proportional to 

r s /r . 

result (Eq. ITHj) without transverse flow. The viscous correction becomes comparable to 
ideal results for pt ~ 1.7 GeV indicating the breakdown of the hydrodynamic description 
of pt spectra for the flow profile considered here. Setting T s /r to 0.1 extends the domain 
of applicability to 2.3 GeV. The analytic Bjorken result (Eq. ITB|) qualitatively explains the 
shape of Fig. El Quantitatively however, the transverse expansion alleviates some of the 
longitudinal shear and pushes the region of applicability hydrodynamics to somewhat larger 
transverse momentum. 

Indeed, viscous effects are implicated in the heavy ion data for px ~ 1.5 GeV. The 
observed elliptic flow deviates from ideal hydrodynamic results for w 1.5 GeV. Further 
for pt ~ 1.5 GeV, the single particle spectra start to deviate strongly from the hydrodynamic 
results (see e.g. j8(). Viscosity provides a simple explanation for the observed breakdown of 
the Pt spectrum in this momentum range. 

Next we examine the effect of viscosity on elliptic flow. In non-central collisions the radial 
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velocity is given a small elliptic component to reproduce the observed elliptic flow 

u r {r , Vs , r, <P) = u ^-(l + u 2 cos(20)) Q{R - r) . (29) 
K 

The functional form of all other hydrodynamic fields is kept the same. Here we simulate the 
STAR 16-24% centrality bin which corresponds to an impact parameter bin (b) ~ 6.8 fm 
jjj. In the model, the radius and lifetime parameters (R Q and r Q ) are scaled downward from 
the central values by the ratio of the r.m.s. radii between b = 6.8 fm and central AuAu 
collisions. This scaling of R Q and r approximates the impact parameter dependence of ideal 
hydrodynamic solutions [8j. The non-central parameters are recorded in Tabled As before, 
once the flow fields are specified, the viscous correction is found by differentiating (^7 a vP). 
The full form of the correction is given in Appendix [X] 

The elliptic flow as a function of transverse momentum V2{pt) is defined by Eq. ^ Ex- 
panding to first order 

/ fdd) * Nm \ f dd>coa(26) d2NW 
V2( P t) = v\ >(p T ) 1 - d2Nm + d2N(0) , (30) 

\ J t p T dp T d<j> / J Pt dpx d(f> 

where {pt) denotes the elliptic flow as a function of px calculated as in Eq. but with 
the ideal distribution ' /A 



Pt ' 



Fig. El shows the elliptic flow for pions. By construction, the ideal curve roughly 
reproduces the experimental elliptic flow at b m 6.8 fm. Taking a more realistic flow profile 
would improve the agreement of the ideal results with data 19||. The effect of viscosity is to 
reduce the elliptic flow. Similar results were recently found [32] by considering a partially 
thermalised expansion. Taken at face value these results suggest that the viscosity is small. 
Indeed, in order to agree with the ideal results up to pt ~ 1.0 GeV we require T s / t < 0.1 . It 
must be mentioned that the results of Fig. El are sensitive to the blast wave parameters. Ideal 
hydrodynamics generates an appropriate set of parameters. Whether a viscous expansion 
(with T s /t = 0.1) can reproduce the observed a elliptic flow remains an open question. 

Finally, I discuss how viscosity effects HBT radii. First, I illustrate the ideal HBT radii 
for the blast wave parametrization in Fig.Efa) The model parameters are again to chosen to 
approximately reproduce the observed radii which are illustrated in Fig.|3fb) for comparison. 
The viscous correction to each radius is again found by substituting / = f + 8f into Eq.l2"4l 
and expanding the numerator and denominator to first order in 5f and calculating the 
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FIG. 3: Elliptic flow V2 as a function of pt for different values of T s /t . The data points are four 



particle cummulant data from the STAR collaboration 



Only statistical errors are shown. The 



difference between the ideal and viscous curves is linearly proportional to T s /t . 




FIG. 4: (a) Ideal blast wave fit to the experimental HBT radii Ro, Rsi and Rl shown in (b) 
as a function of transverse momentum Kt- The solid symbols are from the STAR collaboration 
[~L1 | and the open symbols are from the PHENIX collaboration ^^]. For clarity, the experimental 
points have been slightly shifted horizontally. 
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FIG. 5: (a) Viscous correction SR 2 for Ro, Rs, and Rl relative to ideal blast wave HBT radii 
(R 2 )^ . (b) The HBT radii R , Rs, and R L including the viscous correction. The viscous 
correction is linearly proportional to T s /t q . 

integrals numerically. The resulting viscous corrections are illustrated in Fig.El Several 
observations are immediate. First, as discussed in Sect. IIV[ the viscous corrections in the 
longitudinal directions reduce (z 2 ) and (t 2 ) due to the reduction of longitudinal pressure. 
This reduces the Ro and Rl radii. From a phenomenological point of view the reduction 
of Rl is welcome In full ideal hydrodynamic simulations of heavy ion collisions assuming 



boost invariance in the longitudinal direction 



10, 



33J], Rl is approximately twice too large 



compared to the data. In the blast wave model, viscous corrections to Rl are large. This 
suggests that viscosity is responsible for the shortcomings in these simulations. Comparing 
Fig. 2(b) and Fig. GJb) , it seems that the reduction to Rl is too large. However, it should be 
remembered that the parameters of the blast wave model have been adjusted to reproduce 
the ideal results and therefore viscous corrections make the agreement with data worse. 
Further, because the correction to the longitudinal radius is large the calculation can not be 
considered reliable. For T s /t q w 0.1 the viscous correction to Rl is approximately 30 — 50% 
and the calculation is more reliable. 

Viscous corrections to the transverse variances (x 2 ) and (y 2 ) are small. Consequently, the 
sideward radius receives only a small viscous correction. Viscosity introduces no significant 
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x — t correlation which could influence the ratio of Ro to Rs- In the blast wave model the 
difference between Ro and Rs is due to the contribution (t 2 ). Viscous corrections to (t 2 ) 
are negative and are essentially linearly proportional to this variance. For the particular 
value of T s /t = 1/5 the viscous correction is accidentally correct and makes Ro/Rs ~ 1 as 
illustrated in Fig. E(b). The agreement is accidental but the trend is completely general. 
Viscosity reduces the (t 2 ) and therefore tends to make Rq equal to Rs- This is also welcome 



10 



from a phenomenological point of view. Full ideal hydrodynamic simulations (with 
and without jsJ] the assumption of boost invariance) predict Ro/Rs ~ 1-3 which should be 
compared to ~ 0.9 observed in the RHIC data. 

In spite of these welcome corrections, including viscosity makes some aspects of the 
hydrodynamic description of HBT radii worse. All of the observed radii (denoted generically 
as Rx) scale quite accurately with my = a/ K\ + m 2 as 

Rx oc — L . (31) 



Ideal hydrodynamics readily predicts this scaling (see e.g. |35j, |36j, |37(). Indeed, ex- 
panding Eq. EH for the longitudinal radius of an ideal boost invariant expansion, we obtain 
the Sinyukov-Makhlin formula [35j 

(R 2 L r = rl— . (32) 

ITLt 

Viscous terms immediately break this scaling. Expanding Eq. |2H for the longitudinal 
radius with viscous corrections, we obtain 

^ m+ ^ = ^£-^)- (33) 

Viscous terms break the ideal scaling and this correction grows like ^S- relative to the 
ideal result. This deviation from scaling is not seen in the data. 

There remain several puzzling aspects in the HBT measurements for which viscosity offers 
no explanation. All of the radii are the same order of magnitude and fall with ttit as in 
Eq. El In particular the steep fall with my in the sideward radius was difficult to reproduce 
with the viscous blast wave model described here and in the ideal blast wave model |29(. This 
behavior was predicted based upon a parametrization of ideal hydrodynamics 0, 137^ where 
system cools rapidly during freezeout and where temperature and velocity gradients are much 
larger than the geometric size of the system. It is natural to ask whether these conditions can 
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be dynamically generated from some initial conditions or freezeout dynamics - see [38( for 
efforts in this direction. Large velocity gradients and temperature inhomogeneities should 
increase the relative importance of viscosity. Nevertheless, the success of these models should 
be noted. 



VI. CONCLUSIONS 



In conclusion, I have calculated the first correction to the thermal distribution function of 
an expanding gas due to shear viscosity. The momentum range which is accurately described 
by hydrodynamics is directly related to the shear viscosity and depends upon the particular 
observable. I have estimated this momentum range for single particle spectra, elliptic flow, 
and HBT radii using the boost invariant blast wave model. 

For reasonable values of T s = the viscous correction to the single particle spectrum 

of a blast wave model becomes of order one for pt ~ 1.5 — 2.0 GeV as illustrated Fig. El 

The observed elliptic flow places a constraint on the shear viscosity. Indeed, unless T s /t q 
is less than 0.1, v 2 as a function of px falls well below the ideal curve by px ^ 1.0 GeV. For 
the blast wave model, the viscous corrections to elliptic observables become large before the 
corresponding corrections to the transverse momentum spectra. 

Shear viscosity also plays an important role in the interpretation of the longitudinal 
radius. Indeed, Rl reflects not only the lifetime of the system but also the degree of ther- 
malization in the longitudinal direction. Rl involves the second moment of the thermal 
distribution function in the longitudinal direction where non-equilibrium effects are the 
largest. Consequently, viscous corrections to this radius (approximately 50% for T s /t « 0.2 
and 25% for T s /r « 0.1 .) are large enough that perhaps Rl should be left out of hydrody- 
namic fits to heavy ion data. This does not imply that hydrodynamics must be abandoned. 
On the contrary, while thermodynamics might accurately describe (pt), it certainly does 
not accurately describe (p^ 00 ) unless the viscosity is very small. In addition, viscous cor- 
rections to the ideal longitudinal radius seem to contradict measurements of Rl- Shear 
corrections cause the longitudinal radius to deviate from the scaling clearly seen in the 

nnn n 

data |2ij 120, |3Jj and expected in ideal hydrodynamics [35 1. 

Shear viscosity also reduces the ratio of Ro to R$ by decreasing the emission duration 
(t). Nevertheless, viscosity is not a panacea for the HBT problem. The sideward radius 
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falls precipitously as a function of Kt- This precipitous fall can not be reproduced by 
hydrodynamics at least with a boost invariant expansion 3^|. Viscous corrections to R S ide 
are small and make the sideward radius increase with Kt- 

Many of the conclusions in this work about HBT radii were recently reached "from the 
opposite end" by Gyulassy and Molnar (GM) j^] using kinetic theory. GM, started from 
the Knudsen limit, increased the transport opacity and increased the longitudinal radius. 
Here, I started from the ideal hydrodynamics, increased the viscosity and reduced of the 
longitudinal radius. These authors also emphasized the importance of the y — f] s correlation 
in determining R^. They also found only small viscous corrections to R s and experienced 
similar difficulties in reproducing the steep fall in Kt- 

Clearly performing a full viscous calculation is the next step towards a complete ther- 
modynamic description of the heavy ion reaction. Whether the shear viscosity can be made 
small enough (T s /t < 0.1) in the early stages to reproduce the elliptic flow but still large 
enough (T s /t rs 0.2) in the late stages to reproduce Rl and Ro/Rs remains an open and 
important dynamical question. 
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APPENDIX A: THE VISCOUS TENSOR AND BLAST WAVE MODEL 



To write down the viscous tensor (VaUp) it is most convenient to use Bjorken coordinates: 
r = \Jt 2 — z 2 , rjs = | log (fzf) , r = \Jx 2 + y 2 , and <p = atan(y/x). Note, we denote the 



space-time rapidity with r) s and the viscous coefficient with r\. However, we will drop the "s" 
on raised and lowered space-time indices when confusion can not arise. In this coordinate 
system the metric tensor is 



T 

r 





T 
(\ 








-r 2 





r 



-1 





\ 




-r 2 / 



(Al) 



The only no n- vanishing Christoffel symbols are T^ v = r, = i, = — r, = K 

Without particle number conservation, the hydrodynamic fields are T(r,r] s ,r, <p) and 
u^(r,ri s ,r, 0), where \i = r,T,r] s ,(p. The velocity field satisfies u^u^ = 1 and therefore 
only three components of need to be specified. For boost invariant flow u v = 0. For rota- 
tionally invariant flow = 0. For non-rotationally invariant flow we shall leave = and 
leave the temperature profile rotationally invariant. We assume boost invariance through- 
out. By assumption, the particles freezeout at a proper time r Q with a uniform distribution 
in the transverse plane and a linearly rising flow profile. Thus, the hydrodynamic fields are 
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parameterized as 



T(r o ,^,r,0) = T Q{R -r) (A2a) 

u r (r , n,, r, 0) = u ^-(l + u 2 cos(20)) Q(R - r) (A2b) 

= (A2c) 

u" = (A2d) 



u 



= v 7 ! + (u r ) 2 • (A2e) 



For central collisions w 2 is zero. It is useful to realize that tu v and rw* are the velocities in 
the rj and directions respectively. 

The viscous tensor is constructed with the differential operator V a = A af3 dp, where 
A a P denotes the projector, g a/3 — u a u l3 ) and dp denotes the the covariant derivative, dpu a = 
dpu a +T'^ l3 u fl . With these definitions the viscous tensor is given by rj (V q m^}, where (V a up) = 
V a up+Vpu a — |A a/3 V 7 M 7 . Assuming boost invariance, the spatial components of the viscous 
tensor are given by 

r(vV) = -rd r u* - ^u r - ru r Du* - rvtDvT - -rA^-^d^y^u") (A3a) 

v r 2 1 

r 2 (vV) = -2fyu* - 2 2rVDu* - - r 2 A** -—d^y/^u") (A3b) 



3 

(V r w r ) = -2<9 r w r - 2u r Du r - ^A rr ^=d„(y^u^ (A3c) 

3 v — g 

ii T 2 1 

r 2 ( W> = -2- + --—^(V^u**) (A3d) 
7" 3 V-^f 

(V r M") = <VV') = . (A3e) 
Here y/^g = rr, the expansion scalar is given by 

-MV^g^) = — + — + + d r u r + d T u T , (A4) 



v~g T r 

and the time derivatives in the rest frame Du 11 = u a d a u^ are given by 



Du r = u T d T u r + u r d r u r + u*d^u r - r{ufy (A5) 
rDu^ = w T <9 r (r^) + u r d r {ru 4 ') + d< t> {ru <l> ) + wV . (A6) 

Once the spatial components of the viscous stress energy tensor are known the temporal 
components are determined (numerically) from the relations, (V a u^ up = 0. 
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In these equations the time derivatives, d T u^, d T u r , and d T u T appear. To fix the value of 
these time derivatives it is sufficient to consider the ideal equations of motion. Inclusion of 
viscous terms would lead to previously neglected second order corrections in ^r. The ideal 
equations of motion can be written 

De = -(e+p)V/ (A7) 

Dm" = +— - . (A8) 
e + p 



With these two equations for De and Du r , and the flow profile given in Eqs. IA21 the time 
derivatives can be determined 

d T u* = (A9a) 

C^l) f u T 77'" \ 

d T u r = — + — + d r u r + v 2 d r u r - vd r u r (A9b) 

1 — cr s v l \ t r J 

d rU r = vd T u r . (A9c) 

Here v = u r /u T is the radial velocity and c 2 s = ^ denotes the squared speed of sound. c 2 s 
is very dose to 1 for the pion gas considered and is fonnd by difeentiat.ng the sanation of 
state for a single component massive classical ideal gas. See e.g. |25| for explicit formulas for 
the pressure and energy density. With the necessary time derivatives, the full viscous tensor 
can be found by substituting the flow profile given in Eq. IA2l into Eq. IA3l and differentiating. 
The final formulas are lengthy and are not given. A check of the algebra is provided by the 
trace relation, g^ v T^ s = 0. 

An additional prescription for fixing the time derivative was tried. If the particles are 
freezing out, then the particles are free streaming. Accordingly, we have Du^ = 0. This 
amounts to dropping terms proportional to c 2 s when computing Eq. IA9I This change made 
only a negligible change to final results. This is because the whole effect of the time derivative 
is proportional to c 2 s v 2 which is rather small in practice, c 2 s v 2 ~ 

To finish computing the viscous correction p^p u (V ^Uy) we need to express p^ 
and the integration measure p M cE M in the (r,r) s ,r,<fi) coordinate system. For 
a particle at point (r,T] a ,r,<f)) with four momentum p^ = (E,p x ,p y ,p z ) = 
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(rriT cosh y, pt cos p , pt sin P , my sinh y) we have 

p T = m T cosh(?/ — ?7 S ) (AlOa) 

rp v = rriT sinh(?/ — rj s ) (AlOb) 

p r = pt cos(0 p — (p) (AlOc) 

rp^ = p T sin(0 p - <p) . (AlOd) 

The oriented freezeout volume is dT,^ = (<i£ T , cE r , gE^, dY> v ) = (rdr] s rdr d(j), 0, 0, 0) and the 
integration measure is 

p^dT,^ = rriT cosh(y — r] s ) rdr] s rdr d(p . (All) 

With these formulas there is ample information to compute the viscous correction 
p^p u (V ^u v ) and to perform the necessary Cooper- Frye integrals. 

APPENDIX B: VISCOUS CORRECTIONS TO A BJORKEN EXPANSION 

In this appendix I provide the details leading to the viscous corrections to the spectrum 
and longitudinal radius (Eqs.EI and for a boost invariant expansion without transverse 
flow. The spectrum is given by the Cooper-Frye formula, Eq. El First we compute the 
ideal spectrum. For a boost invariant expansion without a transverse flow u T = 1 and 
vP = u r = = 0. The thermal distribution for an expanding Boltzmann gas is f a (^jr) = 
exp (^ — mT cos ^y~^ a ) ^ Then the Cooper-Frye integral gives the thermal spectrum from an 
expanding cylinder 

^ 1 (Bl, 



d 2 p T dy (2tt) 3 J * V T 

Substituting the integration measure p^dT,^ we have 

^2^(0) i pRo f*r 



1 f u ° f Z7T f°° p. u 

i ^ 7^3 / rdr / d< t> / t drj s rnr cosh(y - T] s ) f (— — ) . (B2) 
d 2 p T dy (2tt) 3 J J J ^ T 

Performing the integral we obtain the ideal thermal spectrum 

d^y= mTT °(2^ 2K ^- (B3) 
Here K\(x) is the modified Bessel function evaluated at x = Now we determine the 
correction spectrum. For a pure boost invariant expansion the non-vanishing components 
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of viscous tensor (V i ii w ) are from Eqs. IA3 



<VV> = (B4a) 

or 

r 2 (vV) = — (B4b) 
3r 

r 2 (VV) = ~ . (B4c) 
or 

Thus the viscous correction 5/ is 

8 T 2 ("r-J P P (V ^ } = 8 T* /o 3r~ Smh * J (B5) 

Note we have substituted / c (l + f Q ) in Eq.[TB]by / as required by the Boltzmann approxi- 
mation. We can then substitute 5f to determine the first viscous correction 

d 2 N^ 1 f 

d^Ty = Wf J ' <B6) 

Substituting the integration measure and performing the integral over the r) s as for the ideal 
case we obtain 

tPprdy (2vr) 3 K Ut\\tJ \tJ yK^x) )) V ; 

Dividing Eq. IB 71 with Eq. IB3I we obtain Eq.lTSl given the text. 

Next we work out the first viscous correction to the longitudinal HBT radius. The 
longitudinal radius is given by Eq.|2Hl Expanding to first order in Sf and using the relation 
z = t sinh r] s we obtain the ideal contribution 



rwfr , = I K ^ US! r ° si nh2 ^ 



and the first viscous correction 



,/V " X f K»dZ, 5 fr 2 sinh 2 Vs 



6R 2 L (K T ) = (R 2 L )W i^L I + j ^ ^^j-o^ ;s (Bg) 



For the kinematics of typical HBT measurements at mid rapidity, we have K 



i-i 



(K T , K r , K* , K") = (y/K% + m 2 ,K T ,Q,0) . The integration measure is K^dH^ = 
vtlt cosh(^ s ) t dr] s rdr d(j) where rriT = a/ Kj, + m 2 . 

First we work out the ideal radius, (R 2 ^)^. Substituting K^dY,^ into the numerator 
and denominator and performing the i nteg rals over the freezeout surface (as in Eq. IB2|) we 



obtain the Herrmann- Bert sch formula |2£ 



^> ,o > - ^'M • (mo) 
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where x = \Jm 2 + K^/T. For large values of x, Eq. IB 101 reduces to the Makhlin-Sinyukov 
formula 35] 

(Rl) (0) = r 2 o — . (Bll) 

A similar calculation gives the viscous correction. Substituting the viscous correction 5f 
(Eq. lB5|l into Eq. lB9[ using the previous results for the spectrum (Eqs. IB31 IB7j) and ideal 
radius (Eq. IBlOj) . and performing the rj s integrals, we obtain Eq. l2T)l quoted in the text 

5Rl T s f6xK 3 (x) x2 l fK 3 (x) \\ (m2) 



r \4 K 2 (x) 8 \K 2 (x) 
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